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We derive from first principles, using non-equilibrium field theory, the quantum 
Boltzmann equations that describe the dynamics of flavor oscillations, collisions, 
and a time-dependent mass matrix in the early universe. Working to leading non- 
trivial order in ratios of relevant time scales, we study in detail a toy model for 
weak scale baryogenesis: two scalar species that mix through a slowly varying time- 
dependent and CP-violating mass matrix, and interact with a thermal bath. This 
model clearly illustrates how the CP asymmetry arises through coherent flavor oscil- 
lations in a non-trivial background. We solve the Boltzmann equations numerically 
for the density matrices, investigating the impact of collisions in various regimes. 
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I. 



INTRODUCTION 



The origin of the baryon asymmetry of the Universe (BAU) is one of the great unsolved 
puzzles in particle and nuclear physics and cosmology. The BAU has been measured through 
studies of (i) Big Bang Nucleosynthesis (BBN), the epoch of light element formation at time 
^BBN ~ 1 minute, and (ii) the cosmic microwave background (CMB), a relic of hydrogen 
recombination at tcMB ~ 10^ years. Characterized by n^/s, the ratio of baryon number 
density to entropy density, the BAU has the value 



at 95% CL. It is a triumph for cosmology that both measurements agree, despite the fact 
that tBBN <^ ^CMB- At the same time, it is a challenge for particle and nuclear physics, as the 
Standard Model (SM) cannot account for the observed BAU [1, 0]. Successful baryogenesis 
requires the concurrence of three necessary conditions j^: violation of baryon number [B); 
violation of C and CP symmetries, where C is charge-conjugation and P is parity; and 
a departure from thermal equilibrium (or a violation of CPT symmetry, where T is time- 
reversal). Therefore, on general grounds a quantitative understanding of baryogenesis mech- 
anisms requires setting up transport equations for quantum systems in out-of-equilibrium 
conditions. More specifically, the genuinely quantum phenomena of particle mixing and fla- 
vor oscillations play an important role in a number of baryogenesis mechanisms, from weak 
scale baryogenesis jsj to leptogenesis 0|, where one needs to follow the evolution of CP 
asymmetries in lepton flavor space js], loj . Mostly motivated by applications to supersym- 
metric electroweak baryogenesis (EWB), in this paper we study the formulation of quantum 
kinetic equations for mixing species from first principle non-equilibrium quantum field the- 
ory. We focus explicitly on mixing scalars in a time-varying mass background, although the 
main concepts and techniques will apply to fermionic systems as well. 

In electroweak baryogenesis (EWB), the baryon asymmetry is generated during the elec- 
troweak phase transition, the era at temperature T ~ 100 GeV when the Higgs field acquires 
a vacuum expectation value (vev). The generation of a net baryon number occurs as fol- 
lows [ij: If the phase transition is first order, bubbles of broken electroweak symmetry 
(($) 7^ 0) nucleate and expand in a background of unbroken symmetry, providing the neces- 
sary departure from equilibrium. CP-violating interactions between the Higgs field(s) and 
other particle species may lead to the production of CP-asymmetries in certain particle 
densities, within the expanding domain wall separating the two phases (bubble wall). These 
CP- asymmetries diffuse ahead of the expanding bubble, and through inelastic scattering 
get partially converted into a CP-asymmetry of left-handed SM doublets {ni ^ 0), which is 
then eventually converted into a baryon asymmetry by weak sphaleron processes. Finally, 
the baryon asymmetry is captured into the broken phase by the expanding bubble, where it 
persists without washout only if the phase transition is strongly first order (($)/T > 1 so 
that electroweak sphalerons are shut off). Within the SM, it turns out that for values of the 
Higgs mass experimentally allowed the phase transition is not first order hence EWB is 
not viable. 

In SM extensions, the viability of this scenario depends on two main considerations. First, 
is the electroweak phase transition strongly first order? Second, is there enough CP- violation 
to generate the observed BAU? 

An answer to the first question follows from a study of the Higgs finite temperature 
effective potential [lo[. The requirement of a strong first-order phase transition in general 



tib/s 



(6.7 - 9.2) X 10-11 BBN \^ 
(8.36 - 9.32) X 10-11 CMB [1 , 2] 



(1) 
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implies strong restrictions on the spectrum. For example, in the Minimal Supersymmetric 
Standard Model (MSSM), this requirement implies a relatively light Higgs boson (rrih < 127 



GeV) and stop (supersymmetric partner of the top quark), with < 125 GeV |lOHl2 



Other models lead to various predictions within an extended Higgs sector |l3Hl6|. The LHC 
should definitely be able to probe these scenarios. 

On the other hand, addressing the second issue remains an open problem that requires 
a quantitative understanding of particle transport at the phase boundary in presence of 
CP violation. In general, one has to write down a network of kinetic equations for the 
relevant species accounting for all relevant processes {CP violation, elastic and inelastic 
scattering, B violation), and solve them under the boundary condition of equilibrium far 
away from the domain wall. Current state-of-the-art treatments utilize the Closed Time Path 
(CTP) formalism [l7l-20| of finite-temperature quantum field theory j21-28||. Yet, within 
this formalism, different groups rely upon different simplifying approximations to compute 
the BAU. Within the MSSM, the most well-studied scenario for EWB, given the same 



underlying parameters, the predicted BAU spans nearly two orders of magnitude [2J, |25|, |28 
Our ability to test the viability of supersymmetric EWB through a combined analysis of 
collider and Electric Dipole Moment experiments j29| depends critically on the resolution of 
this theoretical confusion. 

Within the MSSM, the expanding bubble wall has a thickness L„, ~ 20/T j30[, much 
larger than the typical thermal de Broglie wavelength, the intrinsic length Lint ~ 0{1/T). 
This circumstance enables one to use the so-called gradient expansion in Li^t/L^^ in studying 
the quantum transport dynamics. To leading (first) order in this gradient expansion, the CP 
asymmetries arise through particle mixing [28|. At second order in the gradient expansion, 
one finds a number of additional contributions, including the "semiclassical force" term 
discussed extensively in Refs. j31-35||. Here we confine ourselves to leading order in the 
gradient expansion, in which the dominant non-equilibrium effect is that the bubble wall 
induces space-time dependent particle mixing phenomena for all species coupling to the 
Higgs. Therefore, a remarkably simple picture emerges: physically, the sourcing of CP 
asymmetries is essentially reduced to the phenomenon of coherent flavor oscillation in a 
non-trivial background. On the technical side, mixing particles are to be described by 
density matrices rather than particle number distributions, in order to properly take into 
account the coherent CP-violating oscillations and the collisions that tend to break such 
coherence. The diagonal entries of the density matrix keep track of the populations of the 
individual "flavor states" , while the off-diagonal terms keep track of the coherence between 
the two mixing flavor states. We emphasize that it is only by evolving the full density 
matrices that one can account for the relevant physics and avoid technical issues such as 
the basis-dependence (flavor versus mass basis) of the flnal results. Although the relevance 
of flavor oscillations has been pointed out in Refs. [28|, none of the existing treatments 
in the literature includes correctly the physics of both CP- violating flavor oscillations and 
collisions. 

The present work is the flrst step towards such a complete treatment of the dynamics 
of mixing particles. We study a toy model of mixing scalars ^l,r that incorporates a 
number of salient aspects of the full EWB problem, namely: (i) time-dependent, flavor 
non-diagonal and CP-violating mass matrix, which leads (through flavor oscillations) to the 
generation of a CP-violating asymmetry from a CP-symmetric equilibrium initial state; (ii) 
flavor-dependent collisions in the plasma, described by interactions of $l,_r with a thermal 
background of scalar bosons [A) in equilibrium. Compared to the realistic EWB problem. 
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the main simplification in the toy model is the use of a time-dependent mass matrix, rather 
than one that varies in both space and time. In the latter case, diffusion currents arise and 
tend to enhance the BAU, as they transport the CP asymmetries into the unbroken phase 
where electroweak sphalerons are active 3^. Diffusion currents are not fully incorporated 



in the existing oscillation-based treatment of Ref . [28| . The all- important generalization to 
a space-dependent mass matrix will be the subject of a separate publication, while here we 
focus on the time-dependent case. 

We note that we are not the first to study the problem of flavor mixing in the presence 
of a time-varying mass matrix. In different physical contexts, earlier works 37h41| have 
approached the problem by utilizing the method of Bogoliubov transformations to derive 
equations of motion for particle number operators in vacuum. In the present study, we 
follow a more kinetic theory-oriented approach, as we seek to determine the impact of plasma 
interactions on the evolution of the system (damping of flavor oscillations, equilibration). 



We thus view our results as complementary to those of Refs. |37H41 . 

Our discussion is organized in the following way. In Sec. [TTl we describe our toy model: a 
two-flavor scalar system with a time-dependent mass matrix and interactions with a thermal 
bath, modelled by a third scalar A, assumed to be in equilibrium. Additionally, we provide 
an expansion scheme in the ratios of time scales in the problem. The relevant time scales 
are the collisional mean free time Tcou, the flavor oscillation time scale Tosc, and the wall time 
sea 1^1 r^; we assume that all three time scales are much larger than the inverse frequencies 
of $L,_R, generically denoted Tint- For a slowly varying wall (r^ ^ Tint), this is the physically 
interesting regime. 

In Sec, mil we derive the quantum Boltzmann equation for the two-flavor density matrices 
in our toy model, utilizing the Closed Time Path formalism. Here, our expansion scheme 
plays a key role, allowing for a straightforward generalization of previous one-flavor treat- 
ments. Our results provide a signiflcant improvement over previous treatments of flavored 
Boltzmann equations; we derive the two-flavor Boltzmann equations from flrst principles, 
including collisions, flavor oscillations, and a time-dependent mass matrix in one unifled 
framework. 

In Sec. IIVI we solve the quantum Boltzmann equations numerically for our toy model. 
Previous EWB treatments have relied upon various unproven ansatze for the form of the 
density matrices; in our work, we obtain directly the two-flavor density matrices as functions 
of momentum and time. First, we illustrate of the nature of the (leading-order) CP-violating 
source in EWB, showing explicitly how a CP- asymmetry can arise through coherent flavor 
oscillations in a time- varying background. We show that the resulting CP- asymmetry is 
maximized for Tosc ~ Tw, reminiscent of the "resonant" EWB scenario studied in Refs. 



25j. Next, we explore the impact of collisions with the thermal bath of A bosons. Generally 
speaking, we flnd that collisions lead to decoherence of flavor oscillations and relaxation of 
the density matrices to their equilibrium forms. In the case of flavor- sensitive interactions, all 
CP-asymmetries induced by the wall are ultimately damped away at late times; for flavor- 
blind interactions, CP-asymmetry persists. Lastly, we study how the two-flavor dynamics 
depends on the underlying parameters of the mass matrix. In particular, we flnd that the 
EWB resonance found in previous treatments [21, 22, 2^ does persist somewhat in our more 



We borrow the word "wall" from the EWB problem, even though it doesn't strictly apply to our toy 
model scenario. In this work, "wall" denotes the time region over which the time-dependent mass matrix 
is varying; the "wall time scale" is the time scale over which this variation occurs. 
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exact formalism. 

We provide our conclusions in Sec.|Vl In addition, we provide three appendices. Appendix 
A provides a brief review of the Closed Time Path formalism. Appendix B describes how 
to derive the usual one-flavor Boltzmann equation using the CTP approach. Appendix C is 
an addendum to Sec. IIIIl providing additional technical details related to the derivation of 
the two-flavor Boltzmann equations. 



II. PRELIMINARIES 
A. Toy Model 

We consider a two-flavor scalar system, with fields {^l,^r) = described by the La- 
grangian 

= 9^$^ d"'^ - AP $ + Ant , (2) 

with spacetime coordinate x = =t, x). We couple each flavor scalar field to a real scalar 
A via the interaction 

Ant= (3) 

The A field is assumed to be in thermal equilibrium and we will not consider its evolution 
equations. The role of A in this toy model is simply to provide a thermal bath of scatterers 
for the $ fields. We take the coupling constants 

to be flavor diagonal; indeed, this defines the "flavor-basis" fields $. In addition, we assume 
that the mass matrix M^(t) is flavor non-diagonal and is a function of time: 



m| v(t) e 



-i a{t) 



The off-diagonal elements, assumed to be a function of a time-dependent background field 
(e.g., the Higgs vev], are parametrized in terms of the magnitude v{t) and phase a{t). 
Motivated by EWB [Sy], we assume the following forms 

vit) = ^(l + tanh(t/r„)) (6) 

a{t) = y ( 1 + tanh(t/r„) ) , (7) 

that introduce the external "wall time scale" Tw. A necessary condition for CP-violation is 
d 7^ (see below). Our formalism can be adapted to other functional forms as well, as long 
as the relevant time scale is sufficiently long. 

The important ingredients of this model are as follows: 

• time- dependent, CP- violating mass matrix, which leads to the generation of a CP- 
violating flavor asymmetry from a CP-symmetric initial state; 



• flavor oscillations, due to the flavor non-diagonal mass matrix; 
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• flavor-dependent collisions in the plasma, described here by interactions of ^l,r with 
a thermal background of A bosons. 

It is convenient to transform Eq. (|2]) into the "local mass basis." We diagonalize the mass 
matrix with the time-dependent transformation U(t), such that 



U^{t) M\t)U{t) 
with 



ml{t) 
mi(t) 



1 



cos^(t) 
sin^(t) e^'^W 



sin0(t)e-*'"(*) 
cos^(t) 



tan26'(t) 



2v{t) 



nir — m 



2 ' 
R 



a{t) = a{t) . 



(8) 

(9) 
(10) 



This diagonalization defines the mass basis fields 
the mass basis, is 



[/"f^. The Lagrangian, in 



C{x) = d^(f)^ S'^^ - 0^ - 0t s + 0t s - 0t ^2 + £ 



int ) 



(11) 



where = UW and the dot denotes dt- 

One can study the condition for CP invariance of this model. The action is invariant 



CP 



under the transformation 0j — > rji 0{, where |?7j| = 1, if 
This condition implies that E^^ = S^^ = 0. From the explicit expression 



-e- 



e 



i sin^ e i sin 2^e~*'" 
sin —i sin^ 6 



a 



(12) 



(13) 



we see that CP invariance requires a = 0, or equivalently, d = 0. Once this condition is 
satisfied, the one for z 7^ j is easily satisfied by an appropriate choice of 771,2(0") (now a is a 
constant phase). 

Neglecting the interactions, the free field equations are 



92 + + 2 S 9^0 + + S <\){x) 











(14a) 
(14b) 



The operators in square brackets are Klein-Gordon operators, suitably modified due to the 
time-dependent mass matrix. In this basis, the interaction becomes 



(15) 



with Y{t) = WyU. 
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B. Time scales 



The quantum mechanical evolution of a system out of equilibrium reduces to kinetic 
theory, described by Boltzmann-like equations, in the limit that there exists a hierarchy 
between microscopic and macroscopic time scales. In our derivation of the Boltzmann equa- 
tions below, we utilize a perturbative expansion in the ratios of these time scales. 

The microscopic scale is given by the "intrinsic time" corresponding to the inverse fre- 
quencies of 01^2: 



Tint ~ ^1 2 



- [e + ml,)-''' . (16) 

In the thermal plasma of the early universe (with temperature T), we take Tint ~ T^^. (We 
assume for the zero temperature masses mi^2 ~ T, as typically holds for a subset of the 
superpartners in supersymmetric electroweak baryogenesis; for mi^2 ^ T, these particles 
are not active in the plasma.) 

Typical macroscopic scales in the problem are associated with (i) the time-varying mass 
(with time scale r^), and (ii) coUisional processes that lead to equilibration (with time scale 
Tcoii). In the limit r„, Tcoh ^ Tint, the usual gradient expansion is applicable. In the two- 
flavor case, there arises an additional scale: the oscillation time scale Tqsc = 27rAu;~^, where 
Au = l^i — a;2|. 

In our analysis, we work in the regime in which the following ratios are small parameters: 

Tint Tint Tint /'I '7^ 

^wall = 5 ^coll = 5 ^osc = • V / 

Tw Tcoll Tqsc 

We derive the quantum Boltzmann equations to first order in these parameters, generically 
denoted 0{e). This regime corresponds the following physical picture: 

• The time-dependent background field is slowly varying, such that e„aii <^ 1. 

• The two scalars are nearly degenerate, such that eosc 1- 

• The mean free time is sufficiently long, such that ecoii ^ 1- This assumption is realized 
if the coupling constant y is perturbative, as shown in Sec. IIIII 

Within a more realistic model, our assumption eosc -C 1 is not necessarily satisfied, depending 
on the spectrum. However, for a slowly varying wall, we find that the CP- asymmetry is 
maximized for Tosc ~ t„ (shown in Sec. irv]) : therefore, the regime eosc -C 1 is the physically 
interesting case in which to study these CP-violating fiavor oscillation effects. We speculate 
that when eosc ~ 1, fiavor decoherence occurs rapidly compared to the longer time scales 
Tw, Tcoii; therefore, one can utilize the one-fiavor Boltzmann equation for each of the mass 
eigenstates, thereby neglecting the off-diagonal elements in the mass-basis density matrix. 
However, a complete treatment that is valid in both the small and large eosc regimes and 
that would test this expectation, remains an open problem. 

In principle, in the multi-flavor case, one has to worry about an additional time scale 
beyond those discussed above: Tcoh, the time up to which the quantum mechanical coherence 
among different mass eigenstates persists (and therefore up to which flavor oscillations can 
play a role). One can estimate Tcoh by considering the time up to which wavepackets of 
different mass eigenstates produced in the same process are still overlapping. The key 



physical input here is an estimate of the wavepacket length j42|, which in a thermal bath 
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can be taken as the mean free path of the relevant particles E|. Using this input we find 
that Tcoh ~ ''osc/ecou ~ Tcoii/cosc, which is very long compared to all other time-scales in the 
problem. So in the regime under study, kinematical decoherence effects are negligible. 



III. QUANTUM BOLTZMANN EQUATIONS FROM FIRST PRINCIPLES 

We derive the quantum Boltzmann equations using finite-temperature, non-equilibrium 
quantum field theory, provided by the Closed Time Path (CTP) formalism. In Ref. j43| . 
Calzetta and Hu derived in this way the Boltzmann equation for a single real scalar field with 
collisions El- In this section, we generalize their results in two ways. First, we consider flavor; 
i.e., we study the case of two complex scalar fields, including the effects of mixing between 
flavor and mass eigenstates, quantum mechanical flavor oscillations, and flavor- dependent 
interactions with the background plasma. The physics is similar to that of neutrino oscilla- 
tions in a dense medium 44, Second, we consider a time-dependent CP violating mass 



matrix for the two-flavor scalar system. As discussed in Sec. [H this is a toy model for the 
dynamics of baryogenesis during the electroweak phase transition, wherein one can observe 
some of the essential physics without making an unproven ansatz about the form of the 
distribution functions. 

In this section, we derive the Boltzmann equations pedagogically. A summary of our 
logic is as follows: 

• Derive equations of motion for the CTP Green's functions for the scalar field accurate 
to all orders in e. These equations are (i) the constraint equation, which gives the 
quantum spectrum of microscopic excitations, and (ii) the kinetic equation, which 
gives the macroscopic evolution of the distribution of states in the thermal bath. 

• In order to derive the quantum Boltzmann equations at (leading) 0{e), we truncate 
the kinetic equation at 0{e) and the constraint equation at 0{e^). At these orders, 
the constraint equation identifies excitations in the thermal bath by their tree-level 
dispersion relations, and the kinetic equation describes how their distribution functions 
evolve due to the leading nontrivial effects of collisions, the variation of the "wall" , and 
flavor oscillations. The effects of 0{e) shifts in the dispersion relations add subleading 
corrections to the evolution of distribution functions on top of the most important 
effects we include here. (We deal with the 0(e) constraint equation in Appendix ICl) 

• By integrating the kinetic equation over energy, we obtain the quantum Boltzmann 
equations for the quasi-particle and antiparticle two-flavor density matrices. These 
Boltzmann equations, given at 0{e), describe the evolution of this two- flavor system 
in the presence of a time-dependent and complex mass matrix, accounting for flavor 
mixing, oscillations, and collisions with other particles in the thermal bath. 

For the reader unfamiliar with the CTP formalism, we provide a brief review in Ap- 
pendix |Xl In addition, in Appendix |B] we review the derivation of the one- flavor Boltzmann 



^ We thank Boris Kayaerand Pctr Vogcl for discussions on this point. 

^ The authors of Ref. [43| used a derivation based on the 2 Particle Irreducible effective action, equivalent 
to the approach described below. 



9 



equations, following Ref. 43|, to illustrate the implementation of the above logic in a sim- 
pler case. In this section, we begin immediately with the Green's functions and equations 
of motion in this formalism. 



Closed Time Path Formalism 



In zero temperature equilibrium quantum field theory, only the time-ordered propagator 
is relevant for perturbation theory. At finite temperature and away from equilibrium, in the 
CTP formalism, one requires four Green's functions, corresponding to all possible orderings 
of fields along the closed time path (—00, +00) U (-I-00, —00) (l7l-20|. In addition, to study 



the two-flavor system introduced above, we define a matrix of CTP Green's functions in 
flavor- or mass-basis field space. In the mass basis, we will use indices i,j to label the 
components of the Green's functions, which we define by 

/ GUx,y) -G<Ax,y) \ 
\G>{x,y) -G%[x,y) I 



where 



G>Ax,y) = U{x)(j)]{y)) (19a) 



G<.(x,y) = (0j(y)0,(x)) (19b) 

G^(x,y) = eix' - /) G>^ix,y) + eiy' - x°) G<.(x,y) (19c) 

4(x,y) = e{x' - /) G<.(x,y) + e{y' - x°) G>.(x,y) . (19d) 

These Green's functions satisfy the Schwinger- Dyson equations 



G{x, y) = (x, y) + j <fz j <fw G^^^ (x, z) n(z, w) G{w, y) (20a) 
= G^'^\x,y)+ [ d^z [ d^w Gix,z)Il{z,w)G^^\w,y) , (20b) 



where G^^^ is the free Green's function and IT is the self-energy. Furthermore, Eq. fl2U|) is a 
matrix equation in field space; we have suppressed the field indices. 

Our task now is to recast the Schwinger-Dyson equations into the language of kinetic 
theory. (Here, for completeness, we work to all orders in e; in the next sections, when 
we truncate at 0{e), our results simplify considerably.) First, by virtue of the free field 
equations (fT4|) . the free Green's functions satisfy 

(^a^ + m2(x°) + 2S(x°)9,o + S(x)2 + 9,oS(x)) G^''\x,y) = - 1 5\x - y)I (21a) 
G^'\x,y) (^dl + m'{y')-2dyoJ:{y^) + J:{yY-dyO^{y))=-z6\x-y)I, (21b) 

where / is the identity matrix in both field and CTP component space. Using Eq. (12T]) . we 
act on the Schwinger-Dyson equations Eq. (120!) with the Klein-Gordon operator in Eq. (IT^ : 
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= - i 6\x -y)I-i jd^z n(x, z) G{z, y) (22a) 
G(x, y) ( + m2(/) - 2 S,o S(y°) + S(y)2 - 9,oS(y) ) 

= - i 5\x -y)I-i jd^z G{x, z) I{{z, y) . (22b) 

The Green's functions Eq. (fT9l) contain information both about the spectrum of excita- 
tions of the fields at the microscopic level, and the distribution of states in the thermal 
bath at the macroscopic level. We can decouple the micro- and macroscopic dynamics by, 
first, defining the average and relative coordinates 

X = (X°^t,X) = ^, r^x-y, (23) 

and, then, Fourier transforming with respect to the relative coordinate r. This procedure 
gives the Wigner transform: e.g., 

G{k; X) = j d^r e G(x, y) . (24) 

At zero temperature, the position-space Green's functions only depend on r. At finite tem- 
perature and away from equilibrium, the average coordinate X, describing the macroscopic 
evolution of the system, also plays a role. In this study, we assume isotropy and homogeneity, 
so that quantities depend on t, but are independent of X. 

By taking the Wigner transform of the sum and difference, respectively, of Eqs. f l22p . we 
obtain the constraint equation for G<{k;t): 

2^2 _ ^ ^ G^(A;; t) = e"'* { m^{t) - 2i S(t) + S(t)2, G^}+[ S(t), G^ ] (25) 



+ 1 {n\ } + z {n^ G'} + ^ [n>, g<] +^[g>, n<] 



and the kinetic equation: 



2kodtG^{k;t) = e-'^(^-i[m\t)-2ik^E{t) + E{tf,G^]-i{t{t),G^} (26) 

+ [^^ G^] + [n^ G'^] + i {n>, G<} - 1 {n<, g>}^ , 

where [■,■] and {■,■} denote the commutator and anticommutator. The diamond operator 
O acts on pairs of Wigner transforms according to the definition 

o(.(..,.(..,)^ 1(^11 -1^1^). (.) 
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This differential operator arises due to taking the Wigner transform of the spacetime con- 
volutions appearing in Eq. ( 122|) . In Eqs. ( 125|1 and (126|) we have introduced the combinations 
of CTP Green's functions and self-energies 

= (G* - G*) , n'^ = ^ (n* - n*) , (28) 

where contains information only about the microscopic spectrum, and II'^ gives the 
shift in the dispersion relation for excitations in this spectrum due to interactions. See 
Appendix [C] for further details. So far, the equations of motion Eqs. ( 125|) and ( |26|) are valid 
to all orders in e^aii, Ccou, Cqsc- 

B. Lowest-order solutions of constraint and kinetic equations 

In solving for G^, we utilize a power counting in e. In this section, we work to zeroth 
order in e. First, we summarize the rules for how to determine the order of any term in 
Eqs. ( l25|) and ( l26l) . Recalling that Tint ~ l//^o ~ ^/^i,2{k) and the definitions of e's in 
Eq. f[T7|) . one has that: 

• Each derivative dt acting on U{t) or «^^2(^) carries one power of l/r^. As a conse- 
quence, /cqS scales as k"^ x C(e„aii), S and T? as x C(e^aii)' 

• Each factor of the self-energy n scales as A;QxO(ecoii) (recall that I/tcou ~ li/ko). When 
we evaluate the collision term explicitly below, it will be clear that this expansion is 
equivalent to an expansion in the coupling constant y. 

• Each factor of Am^ = {ml — ml) scales as (cui + U2)'^ x O(eosc)- 

• Each derivative acting on G^{k] t) carries a power of e, up to factors of k^. 

The last rule follows from the kinetic equation (l26ll . truncated at 0{e^). In particular, note 
that the commutator term 

-.[m=W.G^(M)]^-.w(_g,°,_/ff")' 

is C(eosc)- Using the first three rules, the 0{e^) kinetic equation is 

2kodtG^{k;t) = 0; (30) 

therefore, we see that dtG^ is 0{e). Similar arguments at higher order imply that each 
derivative acting on gets an additional power of e. 

We now turn our attention to the constraint equation. At 0{e^), Eq. fl25l) becomes 

(P-m2) Gf^{k;t) = 0. (31) 

At this order we can replace mi, 2 by the average mass m = (mi + 7712) /2, since deviations 
from this replacement are 0{eosc)- In other words, the Green's function G^{k;t) vanishes 
unless the appropriate dispersion relation 

k'^ = ±Uk = ± V|k|2 + m2(t) (32) 
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is satisfiedjl 

A general form of the solution to the constraint equation is then 

Gr;l-'\k;t) = 27r6{e-m') [e{k'){6., + f,,{k,t)) + e{-k')U-k,t)] 

G^^'\k;t) = 2n5{e-m') [^(A;0)/,,(k, t) + + /,,(-k, t)) ] , 

The constraint equation has determined the spectrum of excitations present in the thermal 
bath. We identity fij, fij as the particle and antiparticle density matrices; here, in the free 
field case, they are given by the expectation values and (blbj) of free particle and 

antiparticle mode operators. The evolution of the density matrices is determined by the 
kinetic equation; at 0{e^), Eq. ( 130|) tells us that they remain static. 



C. Kinetic equation at 0{e) 

We now return to the equations of motion Eqs. fl25|) and (1261) for G^ and solve them at 
0{e). We present the solution of the constraint equation at 0(e) in Appendix O Here, 
we are interested in finding the leading nontrivial evolution of distribution functions fij, fij, 
which occurs at 0{e). For this purpose we will need to know the spectrum, or solution of 
the constraint equation, only at C(e°). The solutions to the equations of motion Eqs. (125|) 
and ( 126|) then take the same general form as in Eq. ( l33i) . However, the solution of the kinetic 
equation Eq. ( l26i) at 0{e) will now give the nontrivial evolution in time of the distributions 
fij, fij due to interactions, to the time- varying "wall", and to flavor oscillations. 

Let us now truncate the kinetic equation Eq. (126!) to 0{e). We note immediately that 
we can drop the terms containing and which are 0(e^aii)- Then, we may use 

the information from the previous order that dt acting on Green's functions introduces a 
suppression by at least one power of e. Therefore we may drop all O operators in those terms 
already containing an explicit factor of 0{e), namely, S or H. After these simplifications, 
the terms surviving in the kinetic equation Eq. (l26l) are 

2ko dt G^{k; t) + I [m^{t), G^{k- 1)] + - {m^{t), ^oG^} + 2 A;° [S(t), G^{k] t)] 

1 (34) 
= [li\G^<] + [li^,G^] + ^({H>,G<} - {n<,G>}) . 

There is, however, one remaining term of 0{e^) that can be dropped. In the second term 
on the right-hand side, G^ acts as a source for G^. Since it multiplies a factor of H, which 
is 0{e), we can evaluate G^ at C(e°). Since G^ = G^ — G*, at tree-level we have 

Gj=^Pv(p47;^)5.,, (35) 

which is diagonal. Thus [H^,G''] is proportional to Am^ and so is suppressed additionally 
by Cose- Therefore the term is overall C(ecoiieosc) and can be dropped from the 0(e) kinetic 
equation. 



^ There is a further subtlety if one works in the hmit e-waii, Ccoii ^ 1 limit, but for finite eosc- One finds four 
different modes for the off-diagonal Green's functions: = ±(wi + uj2)/2 and = ±Aa;/2. The latter 
modes become tachyonic with large |k|. However, we find that they decouple from the 0{e) Boltzmann 
equations in the limit eosc <IC 1. A general analysis for tosc ~ C'(l) will have to account for these modes. 
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After suitably rearranging the terms in the kinetic equation, we obtain the master equa- 
tion on which we rely in the remainder of the analysis: 

2kodtG^{k;t) + i [m\t) + iUh{k;t) - 2ik'^j:{t) , G^{k;t)] + ^{m\t),dkoG^{k-t)} = C, 

(36) 

where the collision term C is defined 

Cik;t) = ^({n>,G<} - {ll<,G>})ik;t) . (37) 

On the left-hand side of Eq. (136|) . the term illh shifts the tree-level mass to the interaction- 
corrected value (this is true also in the spectral function, cf. Eq. (1C7I) ). As we will see, 
the S term induces flavor mixing during the time the mass is actively varying, described by 
Eqs. (IMZD- The last term on the left-hand side will disappear in the analysis below since it is 
a total derivative in /cq, over which we will integrate. The collision term C on the right-hand 
side causes equilibration of the distributions /, / through interactions amongst particles, 
antiparticles and A bosons in the thermal bath. 

The kinetic equation Eq. f l36|) is the master equation governing the quantum field theoretic 
evolution of the distribution of states in the thermal bath, expanded consistently to linear 
order in the perturbations of order ewaii, ecoii, Cosc- In the next section, we derive the quantum 
Boltzmann equations for the evolution of the density matrices /, / themselves directly from 
the quantum field theoretic kinetic equation Eq. fl5Bl) . 



D. Prom Kinetic Equation to Kinetic Theory 

In this section, we show how the kinetic equation (15^ has the structure of a quantum 
Boltzmann equation. Using Eq. ( 1331) . we identify the following positive and negative fre- 
quency integrals (which project out the particle and antiparticle modes) as the particle and 
antiparticle density matrices: 

^■(k,t)= / 2A:°G<.(fc,t) (38a) 

dk^ 

/.,(-k,t)= / — (-2fc°)G>.(A:,t) . (38b) 

We obtain the Boltzmann equations by taking the positive and negative frequency integrals 
of Eq. (1361) . Taking the positive frequency integral of the < component of the kinetic 
equation, we have 

/ ^2k'dtG<{k,t)= ^ (-^[m\t) + ^U\G<]-2k'^[J:,G<]+c) . m 

Jo 27T Jo 27T \ ^ J 

Since G'^lk, t) vanishes at the boundaries, the derivative term dG^/dk^ in Eq. fl36l) integrates 
to zero. Using Eqs. ( l38l) . we can express this equation in terms of / and /. In particular, 
there is a trick for evaluating the commutator term: 

f ^ G<] ^ ['-^ (11) G<] ^ /(k.t)| . (40) 



14 



In the first step, we have inserted a factor of unity (since the spectral function, restricted 
to positive frequencies, implies = w); in the second step, we have used Eqs. (!29l [38!) and 
defined 

Since the term Eq. (HOll is already 0{eosc), it suffices to use the zeroth-order constraint 
equation. All in all, the Boltzmann equation for the particle density matrix is, to 0{e), 



dt 



-z[a;,.(t)-zS(t), /(k,t)]+ / ^(C+[^^G<]). (42) 



Similarly, by taking the negative frequency integral of the > component of the kinetic equa- 
tion, we obtain the antiparticle Boltzmann equation 



dt 



I [Mt) + tmJiKt)]- J ^ {C+[U\G>]) . (43) 



In the next section, we evaluate in detail the remaining interaction-induced terms under the 
integrals. 

The simplest physical application of these Boltzmann equations is the case of vacuum 
flavor oscillations. If we neglect both interactions and the S term (induced by the time- 
dependent mixing matrix), the particle Boltzmann equation simplifies to 

^ = -^[u^k,f], (44) 

which is the familiar density matrix formulation of the Schrodinger equation. Given an 
initial condition /(k, 0) at t = 0, the solution to Eq. is 

^/2i(k,0)e^^-* /22(k,0) ) ■ ^^^^ 

To be concrete, let us compute the oscillation probability for — )■ given a pure 
initial state. We define the flavor projection operators Pl = diag(l, 0) and Pr = diag(0, 1). 
In terms of the mixing matrix U, defined in Eq. ([8]), the pure L initial state corresponds to 
the mass-basis density matrix /(k, 0) = W Pl U . The oscillation probability is 

VL^Rit) = Tt[Pr U fit) t/t] = sin^ 26 sin^ , (46) 

as expected. 



E. Interactions and Collisions 

In this section, we evaluate the effects of the interaction 

/:int= -^A2 0ty0, (47) 

in the kinetic equation Eq. fj36|) . This simple interaction, described in Sec. [Tll is a toy model 
for flavor-dependent gauge and Yukawa interactions that are active in the plasma during 
the early universe. 
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FIG. 1. Leading- order self-energy graphs that induce the collision terms in the Boltzmann equations, 
corresponding to (a) coherent forward scattering, and (h) non-forward scattering (cpA (l)A) and 
annihilation (cfxj)^ o AA). 



1. Plasma- corrected mass 



The Uh term in the kinetic equation ( l36l) describes coherent forward scattering, analogous 
to the MSW effect for neutrino flavor mixing j46|. The leading contribution to arises 
from the one-loop diagram shown in Fig. [T](a). The self-energy from this diagram is 



U{x,y) 



--Y{t)6\x-y) 



G\{x,y) 
G\{x,y) 



(48) 



where Ga is the A boson Green's function. Assuming that the A bosons are in thermal 
equilibrium, their Wigner transformed Green's functions are 

G<{p) = (27r) - m\) sign(p°) (49a) 

G>M = (2vr) 5(/ - m\) sign(/) (1 + nB(p°)) . (49b) 



Therefore, we have 



Il\k,t)= --Y{t) 



^ {G>Ap) + G<{p)) 



(50) 



with energy Sp = ^/p'^ + m\. The Ilh term effectively shifts the mass in the oscillation term 



m^{t), G^{k,t) 
In the -C T limit, this shift is given by 

iU''{k,t) = Y{t) 



m\t) + iU''{k,t), G^{k,t) 



d?p 1 



rp2 



(51) 



(52) 



(27r)3 4£pJ 24 

The term in square brackets is the usual zero temperature divergence that must be removed 
through the renormalization of m? . The remaining temperature-dependent term is the one- 
loop thermal mass. This term is flavor diagonal and induces the shift: 



m 



L,R 



m 



L,R 



24 



(53) 



In the limit that iU^ ^ flavor oscillations do not occur since the mass basis coincides 
with the flavor basis. 
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2. Collision term 
Now we evaluate the collision term Eq. fl37|) : 

C^^ (^{U>{k,t),G<{k,t)}-{U<{k,t),G>{k,t)}^ . (54) 

This term encodes the effects of the processes of annihilation (00''' o AA), (non-forward) 
scattering {(pA -H- (pA), and emission and absorption (0j ■H- (pjAA) on the Green's functions 
G^. As we observe below, the latter processes are kinematically forbidden in the small eosc 
regime in which we work, so we will henceforth retain only the annihilation and scattering 
terms. The leading contribution to Eq. fl5^ is the two-loop graph in Fig. [T](b). The resulting 
self-energy is 

U\x, y) = -l Yix'^) G\x, y) Y{y') G\{x, y) G\{x, y) , (55) 
with A =<, >. At linear order in e, its 'Wigner transformation is 

I^Ji^Jw)^^'^^'''^'-''-"-"'^ '''' 

X Y(t)G\k',t)Y(t)G),{p)G\(p') ■ 

Through this self-energy, the collision term's contribution to the particle Boltzmann equation 
(gg) is, using Eq. 




X ( {Yit)G>{k',t)Y{t), G<{k,t)} G>{p)G>Ap') 

-{Y{t)G<{k',t)Y{t), G>{k,t)} G^{p)GW)) . (57) 

We express this in terms of / and / by performing the frequency integrals. Under our ap- 
proximations, there are three integration regions receiving nonzero support from the energy- 
momentum conserving delta function 6'^{k — k' — p — p') that correspond to following three 
physical processes: 

annihilation: 0(fc) (j)\k') ^ A{p) A{p') , (A;°, p'° > 0; k'° < 0) (58a) 

scattering: 0(A;) A{p') o 0(A;') A{p) , (A;°, A;'°, / > 0; < 0) (58b) 

scattering: 0(A;) A{p) o 0(A;') A{p') , (fc°, p'° > 0; / < 0) (58c) 

Other regions of integration correspond to processes that are kinematically forbidden, such 
as (p^k) o 4>{k') A{p) A{p'j^. Only annihilation and scattering are kinematically allowed at 
0{e). 



^ This process is forbidden because, at 0(e), we treat the cj) masses as degenerate within the coUision term, 
consistent with our expansion. At higher order in e, this process does contribute to the colhsion term. 
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Performing the frequency integrals over the region in Eq. fl58al) . we obtain the contribution 
to the colhsion term from annihilation processes, 

-7^C,nn{k,t) (59) 
1 1 fd^k' 1 fd^p 1 fd^p' 1 ^4e4., n 



X {r(t)/(k',t)F(t), /(k,t)} (l + nsM) (l+nB(£pO) 
-{F(t)(l + /(k',t))F(t), (l + /(k,t))} 

Here, the frequencies that appear in the delta function above are all positive; e.g., k'^ = Uk', 
p^ = Ep, etc. We have assumed the A bosons are in equilibrium. 

Similarly, performing frequency integrals over the regions in Eqs. (I58bl) and f l58cp . we 
obtain the contribution to the collision term from scattering, 

°° dk° 

i^Cscat(A;,t) (60) 
Ztt 

---- fA- fp^,^i2.mk-k'+p-p') 

2 2uk J (27r)3 2uk' J (27r)3 2Ep J (27r)3 2ep, ^ ^ ^ ^ ^ ^ 

X ( {Y{t) (l + /(k',t)) Y{t), /(k,t)} UBiSp) (l + nBiSp,)) 

-{Y{t)f{k',t)Y{t), (l + /(k,t))} (l + nBiSp)) nsisp, 
If we take the one-flavor limit of these collision terms, they reduce to the usual semi-classical 



collision terms 43 



We can evaluate these collision terms further. If we assume that the distribution functions 
/ and / are independent of direction of k H, these collision terms take the form 

f°° dk^ 

C[f,f]^ -^Cik,t) (61) 

= -- / dk' Ts{k,k') {l + nB{uJk-ook'))>^ 
Jo '- 

F (1 + f{k', t)) F, f{k, t) } - { Y f{k', t) r, (1 + /(fc, t)) } ) 

+ Ta{k, k') (1 + nsiuJk + ^k')) x 
Y f{k\ t) y, /(fc, t) } - { F (1 + f{k\ t)) F, (1 + f{k, t)) } ) 



with A; = |k| and k' = |k'|. The term with Ta corresponds to the annihilation collision term 
in Eq. (!59l) . while the term with Tg corresponds to the scattering collision term in Eq. (!60ll . 



^ In the case of time-dependent external perturbation considered here, assuming isotropy of the density 
matrices is wcU justified. This assumption will have to be dropped in the realistic case of space-time 
dependent external field. In that case, one has to keep in /(k) the full dependence on |k| and the cosine 
of the angle between k and a unit vector normal to the planar propagating bubble wall. 
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For both annihilation and scattering processes, we have cast the colhsion integral in a form 
that emphasizes the relation between loss and gain terms implied by the principle of detailed 
balance. The anti-particle collision term is simply 



y^^c(fc,t)=c[/,/] 



(62) 



i.e., given by the same expression as Eq. flM]) . but with f ^ f. The rates Ta^g are given by 

J^l pk+k' 



2hQ'K^QkU)k'k j\k-k'\ 

k'T r^+^' 

l2S-K^u)kU)k'k J\k^k'\ 



ds log 
dt I log 



nB(s+)nB(-s^ 



nB(t+] 



to 
T 



where 



So = ^k + ^k' , 
^0 = ^k — ^k' , 



t± ^ ±| + ^Vl-4mi/(t2-t2) 



(63a) 
(63b) 

(64a) 
(64b) 



IV. SOLVING THE QUANTUM BOLTZMANN EQUATIONS 

In the preceding section, we derived the quantum Boltzmann equations for particle and 
antiparticle density matrices {f{k,t) and f{k,t), respectively). The diagonal entries corre- 
spond to the distributions of the individual mass eigenstates, while the off-diagonal elements 
measure the quantum coherence between them. As we have emphasized earlier, it is only 
by evolving the density matrix as a whole that one can properly study the coherent CP- 
violating oscillations (responsible for generating CP asymmetry) and the collisions that tend 
to wash out this coherence. To summarize, our results are: 



df{k,t) 

dt 
dfjk, t) 
dt 



-I [L0k{t) + 5uk{t)-int). f{k,t)]+ g,sC[fJ] 
I [uk{t) + 5uju{t)+int). fik,t)] + g,sC[fJ] . 



(65a) 
(65b) 



The various contributions to these equations are as follows: 

• the flavor oscillation terms [u, f] and [u, /], as discussed in Eq. ( H4l 



the forward scattering terms, giving rise to medium-induced mass terms (analogue to 
the MSW effect), with 



6uJk(t) 



Y{t) f d^p HBiep) 



2oo, J (27r)3 2e, 



(66) 



As previously discussed in Eq. fl53|l . these terms shift the zero-temperature flavor 
diagonal masses: mlj^ — m|^^ -|- ul^rT'^ /2A. Since rriL^R ~ 0{T) and yL,R ~ 0(1), 
we can safely neglect these contributions in the rest of the analysis. 
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• the terms [E, /] and [S, /] arising from the presence of a time-dependent mixing 
matrix U(t); 

• and the colhsion terms C[f, f] and C[/, /], given by Eq. (16T|) . which arise from scat- 
tering and annihilation processes of the mixing scalars ^l,r with a single real degree 
of freedom (A). In Eqs. (!65|) we have introduced the coefficient Qcs, counting the 
number of A-tjpe degrees of freedom present in the thermal bath. In the electroweak 
plasma there are (9(100) degrees of freedom from which the scalars can scatter, so 
Qcs ~ O(IOO) will provide a realistic estimate of the collision rate. In our numerical 
explorations we will vary g^s to dial the relative size of the inverse collision rate Tcou 
versus the other scales in the problem, namely and Tosc- 

In this section, we describe numerical solutions to these equations in a variety of regimes. 
Treating this as a toy model for the dynamics of electroweak baryogenesis, we wish to study 
how the time- dependent mixing matrix (encoded by S) leads to the generation of a CP- 
asymmetry in an initially CP-symmetric plasma, and how this asymmetry evolves in the 
presence of flavor oscillations and flavor- dependent interactions. To illustrate these points, 
we first analyze the collisionless case, paying special attention to the generation of the CP 
asymmetry. We subsequently consider in detail the effect of collisions. 



A. Collisionless flavor oscillations and CP violation 

1. Analogy with spin precession in time- dependent magnetic field 

Neglecting the collision terms, the Boltzmann equations ( 165|) describe the independent 
time evolution of fij{k,t) and fij{k,t), with no coupling between different momentum bins 
or between particles and antiparticles. Without loss of generality, we can parametrize the 
2x2 mass-basis density matrix f(k, t) in terms of four real functions by expanding in the 
identity matrix and Pauli matrices |45|: 

f{k,t) = Ipo{k,t) + a-p{k,t). (67) 

Here, pQ{k,t) is the total occupation number of momentum mode k, and the polarization 
vector p{k) = {px,Py,Pz) describes the density matrix for the "internal" (flavor) degrees 
of freedom. Similarly, we parametrize the anti-particle density matrix fij{k,t) in terms of 
Po{k, t) and p{k, t). 

In the collisionless case, Eqs. (165|) imply that po{k,t) and po{k,t) do not evolve in time, 
while the polarization vectors p{k, t) and p{k, t) obey equations of motion suggestive of spin 



^ Strictly speaking, in this theory, we should then multiply the medium-induced mass correction Eq. 
by 5cff as well. However, in a gauge theory like the electroweak theory or QCD, the leading thermal 
correction to, say, a fermion mass is summed over the number of gauge boson degrees of freedom in the 
theory, whereas the collision terms receive contributions from every particle with which the fermion can 
scatter, such as fermions of other flavors and colors as well as all the gauge bosons. This number of 
degrees of freedom is typically much greater than the number contributing to the leading correction to the 
thermal mass. We mimic this in our simple toy model by the heuristic device of multiplying the collision 
term by geS but not the thermal mass correction. 
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FIG. 2. Left panel: mass eigenvalues mi^2{t) a function of time. Right panel: components of 
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precession: 



dp{k, t) 

Jt 
dp{k, t) 

Jt 



Boik,t) + Bj:{t)j xp{k,t) 
Bo{k,t) - Mt)) xp{k,t) 



(68a) 
(68b) 



with effective time-dependent magnetic field dictated by the form of the mass matrix: 

Bo{k, t) = (o, 0, uji{k, t) - U2{k, t)^ (69a) 
i?E(t) = (2sincr ^ + sin26' cosad, — 2 cos cr ^ + sin 2^ sin acr, 2sin^6'cr) . (69b) 



The particle polarization vector p{k,t) precesses around Bq + Bj^, while the anti-particle 
polarization vector p{k, t) precesses around Bq — in the opposite direction: this describes 
the flavor oscillation dynamics. The usual vacuum oscillation is recovered in the case of t- 
independent mass matrix, in which the effective magnetic field [Bq) points in the ^-direction. 
An initial flavor eigenstate — having px^y{k,t) ^ — precesses around Bq with period Tqsc = 
2n/{uji — UJ2)] see also Eq. (H5|l . 



2. Adiabatic and non-adiabatic regimes 

In the time dependent case, the qualitative behavior of the solution depends on the ratio 
of the oscillation time Tosc = 27r/(a;i — U2) and r^, which controls the time variation of the 



= lO/T vo = ao = it/2 nii = 2.2 T mR = 2T 



TABLE I. Baseline input parameters for the toy model. All dimensionful parameters are expressed 
in units of the temperature T or its inverse. In addition, we take ra^jT <^ 1. 
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effective magnetic field through mi 2, 9, and a. In Fig. [2l we show the time- dependence 
of mi_2(t) (determining -Bo(^) and Tosc) and Bj^if) for our baseline choice of parameters, 
reported in Table [B Let us consider the evolution of a CP invariant, purely L-handed 
initial state, given by {po,p) = {po,P) = (1/2,0,0,1/2) at initial time tin < — Tw- Before 
the external wall turns on (t < — Tw), both Bo and p point along the z axis and there is no 
precession. As the "wall" turns on, the non- vanishing Bj^ tends to push p out of its original 
stationary state, triggering the precession around the time- dependent field Bq + Bj^. In the 
adiabatic regime (r^ ^ Tosc), the polarization vector p effectively tracks the magnetic field 
(with a small precession amplitude that vanishes in the Tosc/tw — ^ limit). On the other 
hand, in the non-adiabatic regime (r^ < Tqsc), when the magnetic field changes on time 
scales comparable to or faster than the oscillation time scale, the polarization vector will 
lag behind the magnetic field and begin precessing with a large amplitude. The precession 
persists at late time {t > r^) around the final constant magnetic field. 

The onset of the non-adiabatic regime is illustrated in the first column of Fig. [31 in which 
we plot the the time-dependence of the particle polarization p(/c, t) with k = 3T (correspond- 
ing to Tosc ^ 35/T, at t = 0). We show four different values of = 40/T, 20/T, 10/T, 5/T 
(from top to bottom), in order of increasing non-adiabaticity, resulting in increasing preces- 
sion amplitudes. Within any realistic model, r„ is fixed and the adiabaticity is controlled 
by k and mi 2. Larger values of k and smaller mass splittings increase Tosc, thus leading to 
increasingly non-adiabatic evolution. 

3. CP violation 

Having identified the basic features of the particle density matrix evolution, we can now 
turn to discuss the effects of CP violation. As discussed in Section [TTl the dynamics are CP 
invariant if and only if a = 0. In this case, with a time-independent phase redefinition of 
$L R, one can set a = 0. Therefore, Bq + B^, and Bq — B-^ lie on the y — z plane; for any 
CP- invariant initial condition (such that p{ti^) = p(tin)), the evolution implies that 

p^{k,t) = p^{k,t) , py{k,t) = -py{k,t) , p^{k,t) = p^{k,t) , (70) 

or, equivalently, f{k,t) = f^{k,t). Geometrically, this means that the angle between p and 
Bq + By. is equal to the angle between p and Bq — By.. On the other hand, if cr 7^ the 
evolution violates CP. Geometrically, the effective fields Bq + B^ and Bq — By. are not 
confined to the y — z plane but develop opposite x components, so that the angle between p 
and Bq + By differs from the angle between p and Bq — By- At late times {t ^ r^, such that 
By — >■ 0), the polarization vectors p and p precess in opposite directions and with different 
angles around Bq (parallel to the z axis); therefore, CP is violated since Eqs. (1701) are not 
satisfied. The above qualitative discussion is confirmed by the explicit numerical solution, 
as illustrated in the second column of Fig. [3|, in which we plot the the time-dependence of 
^{k, t) for k = 3T (with = 40/T, 20/T, 10/T, 5/T, from top to bottom). 

In light of future applications of this formalism to weak-scale baryogenesis, it is highly 
relevant to quantify the CP asymmetry generated in the diagonal elements of the density 
matrix in flavor basis. In particular, we are interested in the L-handed charge density, 
defined by 

nL{k,t) ^Ti [PLU{t) {f{k,t)-f{k,t)) f/t(t)] (71) 
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FIG. 3. Evolution of flavor polarization and CP asymmetry in the absence of collisions, starting 
from a CP-invariant, pure L-handed initial state. Left column: time dependence of the particle 
polarization vector components px (red), py (green), and pz (blue), for different values of = 
40/r, 20/r, lO/T, 5/T (from top to bottom). Middle column: time dependence of the anti-particle 
polarization vector components Px (red), Py (green), and Pz (blue), for different values of t„ = 
40/T, 20/r, lO/T, 5/T (from top to bottom). Right column: time dependence of the flavor diagonal 
CP asymmetry nL{k,t) for different values of Tw = 40/T, 20/T, 10/T, S/T (from top to bottom). 
In all cases k = 3T and all other input parameters ( except Tw) are as in Table corresponding to 
rose - 35/T (at t = Q). 
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FIG. 4. Maximum value of CP asymmetry \nL{k,t)\ for k = 3T as a function o/Tosc(t = 0)/tw. 
Except for that is varied, all other input parameters are as in Table Q- 

where Pl = diag(l,0). We plot riL^k = 3T,t) in the third column of Fig. [3l as usual for 
different values of r„ = 40/T, 20/T, 10/T, 5/T (from top to bottom). The flavor-diagonal 
CP asymmetry is small in the quasi-adiabatic regime (top panel) and increases with Tosc/tw 
(the evolution becomes more non-adiabatic). 

However, as one increases the non-adiabaticity, a maximal CP asymmetry is reached. 
This is illustrated in Fig. E] by the fact that the CP asymmetry is larger for r„ = 10/T 
than for r„ = 5/T. As a measure of this effect, we plot in Fig. IHthe maximum of \nL{k = 
3T, t)| versus the ratio rosc(^ = 0)/t^. The CP asymmetry vanishes in both the adiabatic 
(tosc/t"w ^ 1) and extreme non-adiabatic (tosc/tw ^ 1) limits. A maximal CP-asymmetry 
occurs for Tosc/tw ~ 4. This result is independent of the particular choice of momentum bin 
k = 3T. 

The behavior of these limiting cases can be simply understood on a qualitative level. In 
the adiabatic limit (tosc/tw ^ 1), both p and p track their effective magnetic fields Bq ± By.. 
Since at late times these effective magnetic fields coincide (-Bs = 0), the resulting p and 
p are equal and no CP asymmetry is present. On the other hand, in the extreme non- 
adiabatic regime (tosc/tw ^1), the time behavior of the effective magnetic field is singular: 
-Bs(i) ~ S{t) bj]. This produces a common shift 6e x p[0) in both p and p at t = 0. After 
this modified initial condition, p and p precess around the constant Bq in exactly opposite 
directions, so again no CP asymmetry can arise. El So for the system to feel the CP-violating 
effects, Ps has to vary not too fast compared to the oscillation scale. 

In summary, we have shown how a CP-asymmetry can arise in an intially CP-symmetric 
plasma in equilibrium through time-dependent, CP- violating flavor oscillations. This pro- 
vides a conceptually simple picture for the EWB mechanism. We have also identifled and 
discussed the regime in which the CP asymmetry is largest, namely Tosc ~ Tw- We discuss 
next how interactions between the mixing scalars and the thermal bath affect this picture. 



Strictly speaking, a time-independent phase redefinition has to be performed to satisfy conditions (j70p . 
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B. Full Boltzmann equations 

1. General considerations 

The primary technical challenge in solving the full Boltzmann equations is that they 
are coupled integro-differential equations. We surmount this obstacle by discretizing / and 
/ into N momentum bins, each of width A. Whereas Eqs. ( !65l) are 8 coupled integro- 
differential equations (each / and / are composed of 4 real functions), discretization converts 
the Boltzmann equations into 8 coupled ordinary differential equations, which can be 
handled numerically. We have explored several choices of binning, keeping in mind two 
criteria: (i) the need to cover the entire thermal spectrum up to sufficiently large momentum 
(we have included up to /c ~ 8T); (ii) the need to ensure thermalization, which constrains 
the bin width to be not much larger than the width of the scattering and annihilation kernels 
^s,a{k, k') (see Eqs. flB5|) ). For our baseline numerical analysis we have adopted the following 
choice: 

fik,t) ^ fnit) = fikn,t) n = l,N = 8; kn=(n-^^A; A = 1.2T , (72) 

and we have verified that the main results are stable at the percent level against changes of 
binning that satisfy criteria (i) and (ii) above. 

Let us now turn to the physics. At late negative times t <C — t„, we assume the system 
to be in thermal equilibrium with zero total charges. Hence, our initial condition is 

fij{k,tin) = fij{k,tin) = j , fii{k,tin) = fii{k,t;^) = nB{uJi{k)) . (73) 

As shown in the previous section, in the non-adiabatic regime the turning on of the wall 
kicks the system out of equilibrium and induces flavor oscillations. The inclusion of C[/, /] 
in the kinetic equations (l65l) has two important effects on the evolution of the system: 



Collisions tend to destroy the quantum coherence characterizing the evolution of linear 
combinations of different mass (or flavor) eigenstates. This is highly relevant to our 
problem, as coherence is essential for flavor oscillations and hence CP violation to 
emerge. Even iso-energetic collisions would stop the coherent development of the 
state, as long as the two mixing particles have different couplings to the thermal bath 



(Vl 7^ Vr in our model) [4J, |45|]. So on general grounds we can conclude that collisions 
will push the off-diagonal elements of the density matrices f{k,t) and f{k,t) (hence 
Px,y{k,t) and px,y{k,t)) to zero. 

More generally, collisions allow for energy exchange between the scatterers and the 
thermal bath, thus pushing the system towards its thermal equilibrium state. As 
a consequence, the diagonal entries of the density matrices tend to the appropriate 
thermal form consistent with quantum statistics and the conserved charges in the 
system. 



2. Reaching equilibrium 



In the toy model under study, the form of the late-time equilibrium state depends on the 
interaction coupling constants yL,R- By analyzing the collision term, one can easily verify 
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that the general form of the equihbrium distributions will be 

fn{k) = nB{u)k - /ii) fii{k) = nB{u)k + l^i) (74a) 

f22{k) = UBiuJk - 1^2) f22{k) = nsiuJk + ^^2) , (74b) 

with vanishing off-diagonals. In the case of flavor-sensitive interactions (^/x, 7^ ur), one also 
has the condition ^2 = f^i, since there is only one conserved charge in the system, the 
total charge Qi + Q2- On the other hand, flavor-blind interactions {i/l = Vr) can lead 
to equilibrium with fii 7^ fi2, since at late times Qi and Q2 are separately conserved. In 
addition, if we begin with a CP-symmetric initial state (with Q1+Q2 = 0), then at late times 
we must have /xi = — /i2- Therefore, we expect our numerical solutions to reach the form 
given in Eq. ( 174|) . with (i) yUi = /i2 = for flavor-sensitive interactions, or (ii) /ii = — /i2 7^ 
for flavor-blind interactions. 

In the flavor-blind case {yi = Vr), we can compute the individual charge generated at 
late times by the discretized momentum integral 

N ,2 

Qiit) = Aj2^{fuik,t)-fn{k,t)) , (75) 
1=1 

and, using Eq. (1741) . the corresponding chemical potential (working to linear order in ^i/T): 

-1 

(76) 

We have explicitly verified that our numerical solutions converge at late time to these ex- 
pected equilibrium results, turning the above relations into a nontrivial check of our numer- 
ical codes. 

3. The effect of collisions: TcoW ^ and TcoU < Tju 

In this section, we study quantitatively the impact of interactions with the thermal bath 
through a set of simulations whose results are shown in Figs. |5] and [61 Our key results are: 

• Collisions (both flavor-blind and flavor-sensitive) lead to damping of flavor oscillations, 
such that px^y, px,y 0. 

• For flavor- sensitive collisions, all charge induced by the wall is damped away at late 
times {pz, — )■ 0). 

• For flavor-blind collisions, charge induced by the wall is not damped away at late times 
{Pz 7^ 7^ 0). This charge can be interpreted as a chemical potential. 

• Fast collisions (tcou < r„) suppress the generation of CP asymmetry. 

Fig. E] illustrates the impact of flavor-blind versus flavor-sensitive interactions; it is or- 
ganized as follows: In the left column of Fig. [5] we plot the time evolution of the particle 
polarization vector p{k = 3T,t), while in the right column we show the evolution of the 
flavor-diagonal CP asymmetry nL{k = 3T,t). We take different choices for yL,R, illustrative 
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FIG. 5. Evolution of flavor polarizations with no collisions, flavor-blind collisions, and flavor- 
sensitive collisions. We plot the components px{k,t) (red), Py{k,t) (green), and pz{k,t) (blue) 
[left panels] and nL{k,t) [right panels], for k = 3T, as a function of time in several regimes: no 
collisions (1st row); flavor-blind collisions, Ul = UR = deff = 200 (2nd row); flavor-sensitive 
collisions with = 1, Ur = 0.8, Qeff = 200 (3rd row); flavor- sensitive collisions with ul = 1, 
Ur = 0.5, g(,ff = 200 (4th row). In all cases we use equilibrium initial conditions at t\^ = — IS/T. 
All other input parameters are as in Table\^ The dotted line in the 2nd row represents the nonzero 
density surviving at late times corresponding to a non-vanishing chemical potential ^l- In the 
3rd and 4th rows, the late-time density goes to zero since yi yR- See text for additional details. 
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of several regimes: (i) top row: collisionless case Ul = Ur = 0; (ii) second row: flavor-blind 
case Vl = Ur = 1; (iii) third row: flavor- sensitive case = I, i/r = 0.8, illustrative of 
a squark-driven EWB scenario, where the two flavor eigenstates share a common strong 
interaction and are distinguished only by electroweak interactions; and (iv) fourth row: 
flavor- sensitive case yi = 1, Vr = 0.5, illustrative of the chargino-like scenario in the MSSM, 
in which the two flavor eigenstates have interaction strengths that might differ by a factor 
of 0(1). We have adopted gcs = 200, corresponding to an inverse collision rate Tcoii consid- 
erably larger than than the characteristic time scale of the external field = 10/T. We 
evolve the density matrix from equilibrium initial conditions at tin = — 15/T, using input 
parameters as in Table [B 

The panels in Fig. [5] clearly illustrate the impact of collisions on flavor oscillation dy- 
namics: in all cases the expected late-time thermalization is reached. In the case y^ = yR 
(second row), while the coherences decay {px,y{k,t) 0) at late times, we find pz{k,t) — )■ 
{riBioJk — fJ'i) — nB{u!k + yUi))/2, with /ii given by Eq. fl75]) . Similarly, the late time behav- 
ior of nL{k,t) = fLL{k,t) — fLL{k,t) is consistent with the equilibrium form for the flavor 
diagonal density fiLik^t) = nuioJk — /^l), with effective L-handed chemical potential given 
by /iL = (cos^ 6^ — sin^ 6*) /ii □. The equilibrium value of n^^ corresponding to this /iL is 
shown as the dotted line in the second row of Fig. [5l Introducing flavor-sensitive interac- 
tions [yi 7^ yR, third and fourth rows) induces no big differences at early times (t ~ Tw for 
QeK = 200), while it modifies the late time behavior as now pz{k,t) — )■ and riLik.t) — )■ 0. 
By comparing the behavior in the third and fourth rows of Fig. [5l one can also observe that 
as \yL — yR\ grows, the chemical equilibrium (/ii = /x^. = 0) is reached faster, in accord with 
intuitive expectations. 

So far we have used values of the input parameters so that Tcou ^ Tw,Tosc- In this case 
the dynamics divides into three regimes, (i) Early times t = 0(rw) ~ 10/T. This is the time 
scale over which the CP-asymmetry is generated through flavor oscillations, and to a first 
approximation one can neglect the effect of collisions, (ii) Intermediate times t ~ 100/T. 
This is the time scale when collisions damp away flavor oscillations. Over these times, it 
appears irrelevant if the interactions are flavor-dependent or flavor-blind, (iii) Late times 
t ~ 1000/T. This the time scale over which chemical equilibrium is reached. Here it matters 
crucially if the interactions are (nearly) flavor-blind or flavor- dependent. 

If the underlying parameters are such that Tcou < r^, Tosc, one expects a strong damping 
of flavor oscillations and a strong suppression of the CP asymmetry. In this case the role 
of collisions is not merely to relax to equilibrium the asymmetry generated by turning on 
the wall. Collisions are now so frequent that they prevent the system from going sufficiently 
out of equilibrium; oscillations cannot play a significant role in generating a. CP asymmetry. 
We illustrate this effect quantitatively in Fig. [6l in which we change Tcou by dialing the 
parameter (yfes representing the effective number of degrees of freedom present in the thermal 
bath. Again, we adopt equilibrium initial conditions at tin = ~15/T, use input parameters 
as in Table [H and adopt yi = ^,yR = 0.5. In the left column of Fig. E] we plot the time 
evolution of the particle polarization vector p(/c = 3T, t), while on the right column we show 
the evolution of the flavor- diagonal CP asymmetry nL{k = 3T, t). From top to bottom we 
have gcs = 200, gcs = 1000, Qcs = 2000, respectively. 



^ At late time /ll can be expressed in terms of /n and /22 via /ll = cos^ fn + sin^ 9 f22- Using these 
relations, and expanding to linear order in /i, one can relate to /ii as described. 
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FIG. 6. Dependence of solutions on the relative size ofr^ and TcoU- Keeping r„ fixed, we control t^oW 
by dialing the effective number of degrees of freedom g^s in thermal bath. We plot the polarization 
vectors px{k,t) (red), Py{k,t) (green), and pz{k,t) (blue) [left panels] and nL{k,t) [right panels], 
for k = 3T, as a function of time with yi = 1, y^ = 0.5 and all other input parameters are as in 
Table\^ for different values of g^s- gcs = 200 (1st row); g^s = 1000 (2nd row); g^s = 2000 (3rd 
row). In all cases we use equilibrium initial conditions at t[^ = —15/T. See text for additional 
details. 



4- Is there a resonance? 



In the realistic MSSM case, the baryon asymmetry depends on the amount of L-handed 
charge asymmetry that survives in the unbroken electroweak phase after being generated 
within the bubble wall. While the toy model, with only a time- dependent mass matrix, 
does not allow us to obtain a realistic estimate of the space-dependent L-handed charge 
distribution, it does allow us to study how the generated L-handed density depends on the 
underlying parameters of the mass matrix. In particular, our toy model can reveal whether 



the so-called "resonant" enhancement found in Refs. mmm for rriL ran survives within 
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FIG. 7. Dependence of the late time L-handed chemical potential fi^ on the mass parameter mi, 
with mji/T = 2 and all other input parameters as in Table\^ with flavor-blind interaction yi = 
VR = 1- 



the full particle mixing treatment. To address this question, we need a measure of the 
total charge asymmetry generated. For this purpose, we find it most convenient to work 
in the flavor-blind limit ijl = Ur and evaluate the effective L-handed chemical potential 
fJ'L = /ii(cos^^ — sin^ ^) that emerges in the late-time solution, starting from equilibrium 
initial conditions at some early initial time tin < —Tw In Fig. [7]we plot the behavior of ^l/T 
versus mi,/T, with all other input parameters fixed to the values of Table |T] (in particular 
m^/T = 2). The behavior shown in Fig. [7] results from the competition of several effects: 

• For rriL = m^, there is no generation of CP asymmetry at all, because the equilibrium 
initial condition with equal masses implies that p{k,tirt) = (0,0,0), so there is no 
effective source (the commutator terms [S, /] and [S, /] vanish). As long as tul 7^ 
rriji, the system gets pulled out of equilibrium when the wall turns on. However, for 
mi — )■ the initial condition Pzik^U^) — ?■ 0, thus suppressing the final asymmetry. 

• On the other hand, for given momentum as mi — )■ the ratio Tosc/Tw increases 
and so one enters the non-adiabatic regime. This tends to increase the CP asymmetry, 
with maximal CP- violating effects obtained for Tosc/tw ~ 4 (see Fig. Hj). 

• The chemical potential /ii is enhanced when the mixing angle is maximal {9 = 7r/4) 
and suppressed for 6* = 0. On the other hand, the resulting /i^ is suppressed when 
6* — 7- 7r/4, due to the cos 26* prefactor. 

It will be intriguing to study the extent to which this double-peaked structure will persist 
beyond leading non-trivial order in ecou (see, e.g., [2I, [2^) or when the bubble wall evolves 
not only in time but also in space. 

In summary, in this section we have studied the effect of collisions on the generation and 
evolution of flavor-diagonal CP asymmetries. In preparation for more realistic studies in 
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supersymmetric scenarios, we have explored different regimes, by varying: (i) the degree to 
which interactions distinguish the two mixing ffavors {yi versus yji); (ii) the ratio t^/tcoh of 
the wall time scale over the typical collision time. In the supersymmetric EWB scenarios, 
several mixing species can potentially contribute to the baryon asymmetry (squarks, staus, 
charginos, neutralinos, or Higgs scalars) and clearly their interactions fall into different 
regimes. While here we do not attempt to make predictions for the realistic MSSM case, we 
can draw two general conclusions: (i) a larger asymmetry is generated as long as Tw/tcou < 1; 
(ii) the generated asymmetry persists longer if yi/VR is close to unity (the case of roughly 
flavor insensitive interactions). 

V. CONCLUSIONS 

In the present work we have presented the first step towards a consistent analysis of the 
transport of mixing particles in a CP- violating external background. While our ultimate 
goal is to apply the insight and techniques to the problem of weak scale baryogenesis, here 
we focused on a simplified model of mixing scalars $l,_r with time-dependent mass matrix, 
capturing the central ingredients of the full baryogenesis problem. We derived quantum 
kinetic equations for mixing scalars from first principles in non-equilibrium quantum field 
theory in the regime in which the oscillation time scale Tqsc is comparable to or longer than 
the external "wall" time-scale r^. This is the physically interesting regime of non-adiabatic 
evolution, in which the non-equilibrium and CP- violating effects are largest. 

Our analysis provides a novel, simple physical picture for the generation of flavor diagonal 
CP asymmetries in the weak scale baryogenesis scenario: starting from a CP-conserving 
equilibrium initial state, a CP asymmetry arises through coherent flavor oscillations. We 
have shown that, fixing the underlying CP- violating phases, the CP asymmetry is essentially 
determined by the ratio Tosc/tw, reaching a maximum for Tosc/tw ~ 4. 

We have also studied in detail the effect of collisions, by solving the coupled kinetic 
equations for particle and anti-particle density matrices in a number of different regimes. 
As expected, collisions lead to decoherence and hence tend to suppress the CP- violating 
asymmetry. From our analysis, two general lessons can be drawn: (i) a larger asymmetry 
is generated as long as t^/tco\\ < 1, i.e. for mean free path longer than the time scale over 
which the wall turns on; (ii) the generated asymmetry persists longer if collisions are nearly 
flavor insensitive (that is, ii yi/yn is close to unity in the toy model). 

To make contact with the existing baryogenesis phenomenology, we have studied the 
dependence of the CP asymmetry on the underlying mass matrix parameters, in particular 
on the difference of the flavor-diagonal mass entries rriL and m^j. Working to linear order in 
the (small) ratios of physical timescales, we found that while the non-adiabaticity condition 
tends to maximize the asymmetry for ~ rriji, in the case of exact degeneracy = 
the asymmetry vanishes. So some sort of "resonant behavior" persists but with a double 
peak shifted above and below the rriL = mn point. 

As already emphasized, the results presented here represent a first step towards a realistic 
treatment of transport phenomena involving mixing particles at the weak phase transition. 
A number of outstanding issues have to be clarified before more realistic BAU calculations 
can be performed. In particular, work is in progress to: 

• Include a spacetime-dependent mass matrix. In this case the main novel ingredient is 
the generation of spatial CP- violating currents. We plan to study them numerically in 
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the toy model and compare to the diffusion approximation. A proper understanding 
of spatial currents (whether or not in the diffusion limit) is crucial for the application 
to baryogenesis, as currents are known to enhance the BAU 361] by transporting CP 
asymmetries in the unbroken phase where sphalerons are active. 

Formulate kinetic equations including effects beyond linear order in e. 



Include in this formalism inelastic (particle-number changing) interactions |47 



Extend the formalism to fermions and extend the network of kinetic equations to 
realistic cases in the MSSM and other extensions of the SM 48-52 . 
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Appendix A: Review: Closed Time Path Formalism 

In this section, we provide a brief introduction to the basics of the Closed Time Path 
(CTP) formalism. 

Let us begin by considering a single complex scalar field governed by the Lagrangian 

C = d^ip*d^ip - m^ip*ip + £int . (Al) 

At zero temperature, perturbation theory is essentially the study of time-ordered propaga- 
tors, such as 

G\x,y) = (r{M^)^liy)}) , (A2) 

where iph is the Heisenberg-picture field. The key difference when moving to finite temper- 
ature is that the expectation value in Eq. (]A2p is calculated not in the vacuum but in an 
ensemble of states defined by a density matrix 

P = ^Wn\nh){nh\ , (A3) 

n 
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where the time-independent Heisenberg-picture states \nh) each have weight Wn- 

Now, let us move to the interaction picture. First, the interaction-picture states \n{t)) 
are functions of time; we define the interaction states at time t = — oo to coincide with the 
Heisenberg states: |n_) = |n(— oo)) = \nh). The density matrix is p = ^„ w„|n_)(n_|. 
Second, the interaction fields ip are related to their Heisenberg counterparts by the time- 
evolution operator U 

iph{x) = U (xo, -oo)"^ Lp{x) U -oo) . 



The operator U obeys the usual relations: 

and ^ 

ty{ti,t2) = r jexp (^i I dz^ ! d^z C-, 

With these relations, Eq. ( 1A2I) becomes 



--intl 



T|exp i j 



d'^zCintiz) 



(A4) 

(A5) 
(A6) 

(AT) 



d zC\r,dz 



where J d'^z = J^^dz^ J d^z. Reading from right to left, Eq. (1A7P corresponds to starting 
with the "in"-state then time-evolving from — oo to -f-oo, acting with the field operators 
at times and along the way, and lastly time-evolving from -|-oo back to — oo, returning 
to the "in"-state. This time-contour, denoted C, is the "closed time path"; it is closed in the 
sense that the contour begins and ends at t = — oo, connecting "in" -states with "in" -states. 
Eq. (1A7P can then be succinctly written as 



G\x,y) = (V^ip4x)ipliy)exp 
V \ (^+(x)(^^(y)exp 



i I d'^zCiritiz) 



I j d'z (^Cj+\z)-Cj~\z)) I 



(A8) 



where V means path-ordering of fields along C. In the second line, we have reverted to a 
time integral from — oo to -|-oo, splitting Cint into two pieces corresponding to the forward 
and backward branches along the closed contour C. The notation ip±{x) and £int'-''^''(^) — 
itself a function of ^p±{x) — denotes whether x^ is on the time- increasing (-I-), or time- 
decreasing (-) branch of C. The path-ordering prescription is to time-order the (-I-) fields, 
to anti-time-order (T^) the (-) fields, and lastly to put all the (-) fields to the left of the 
(+) fields. 

A perturbative evaluation of y) proceeds similarly to zero-temperature field theory. 
Wick's theorem applies as usual, but with "P-ordering instead of T-ordering. Therefore, we 
must consider not one but four different propagators, corresponding to all possible path- 
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ordering of x° and ?/° in {ip[x)ip\y)): 

G>(x, = ( P [ip-{x) ^\{y)]) = { ^{x) v^y) > (A9a) 

G<{x,y) ^ (^v[ip4x)vUy)} ) = {v\y)v{x)) (A9b) 

G\x,y) ^ ^v{Mx)vi{y)]) = {r{y^{x)v^{y)}) (A9c) 

G*>,y) ^ (p{^_(x)^L(y)}) = (7^{^(x)^t(^)|> (Agd) 

These Green's functions are the free or full propagators for fields in the interaction- or 
Heisenberg-pictures, respectively. In particular, we see from Eq. (lASP that a perturbative 
expansion of G^{x,y) will necessarily involve contracting (+) and (-) fields together; these 
additional propagators are essential. The OTP propagators can be assembled into the matrix 



G\x,y) -G<{x,y) 
G>{x,y) -G\x,y) 



] ■ (AlO) 



In this form, the propagators satisfy the relations 

G(r.,y)^ g<«. + (...) St.. u.) .)) (Alia) 

= G(o)(x,i/) + j (fw j d'^z {g{x,z) Il{z,w) G^^\w,y)^ . (Allb) 

These equations are the CTP version of the Schwinger-Dyson equations, where now both the 
propagator and the self-energy 11 (defined by Cint) are 2x2 matrices. The free propagator 
G''^°^(x, y) satisfies 

[dl+ml) &^\x,y)={dl + ml) G^'\x,y) = -z 6\x - y)I , (A12) 

where / denotes the 2x2 identity matrix in CTP propagator space. 

Taking the Wigner transform Eq. (12^ of the Green's functions in Eq. ( lAlip . and then 
taking the sum and difference of the two equations, we obtain the constraint and kinetic 
equations for the CTP propagators. 



^-k' + rn'\G{k;X) = -zr ^) (Al3a) 



^e-o ( {U{k- X)}{G{k- X)} + {G{k- X)}{U{k- X)} 



2k-dxG{k;X) =6-'"^ (^{U{k;X)}{G{k;X)} - {G{k;X)}{U{k;X)}^ , (A13b) 

where the O differential operator was defined in Eq. (1271) . 

The various CTP Green's functions contain information both on the spectrum of exci- 
tations of the if fields in the Lagrangian Eq. flAip and on the distribution of states in the 
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thermal bath. We can isolate Green's functions depending only on the spectrum by forming 
the retarded and advanced propagators, 

G^{x,y) = (G* - G<)ix,y) = - y')iG> - G<){x,y) (Al4a) 

G^{x, y) = (G* - G>){x, y) = 9{y' - x°)(G< - G>){x, y) , (Al4b) 

whose poles in momentum space give the spectrum of ip excitations. We also use the linear 
combinations, 

G\x,y) = ^{G' - G')ix,y) = ^e{x' - y°)(G> - G<)(x,y) (Al5a) 

p(x, y) = {G> - G<)(x, y) = (G« - y) . (A15b) 

The latter, p{x,y), is also called the spectral function. Note in addition the relations 

= G\x,y) ± ^(G> - G<)ix,y) . (A16) 

We will also define analogous components of the OTP self-energy 11. In particular, if we 
define 2ir = IP' — U^, we obtain 

n^'^ = n'^±zr. (Air) 

The n'* component of the self-energy induces a correction to the mass in the dispersion 
relations for ip excitations, while F induces a nonzero width. 



Appendix B: One Flavor Boltzmann Equation 

In this appendix we show how the OTP equations of motion above may be used to derive 
the Boltzmann equations for the simple case of a single scalar field with a constant mass. 



following the logic of Ref. (43|. We consider the one- flavor version of the interaction Eq. 



between if and a scalar field A, 

We will work from the Schwinger-Dyson equations Eq. flAlip and the constraint and ki- 
netic equations Eqs. flA13aP and (IA13bp . solving them order-by-order in an expansion in y, 



and therefore, ecou- We will derive equations for the spectrum of excitations through the 
constraint equation and equations for the retarded and advanced propagators, and for the 
distribution functions from the kinetic equations for G^. At 0{y), we will obtain a correc- 
tion to the spectrum of if excitations, and at O^y"^), we will be able to derive Boltzmann 
equations for the distribution functions from the quantum kinetic equation with a collision 
term driving the distributions to equilibrium. 



1. Tree- level: Generic Form 



At tree level, the self-energies 11 = 0, so the constraint and kinetic equations Eqs. (IA13ap 
and ( 1A13b|) are just 



,2 



r + J G{k- X) = -I 
k-dxG{k;X) = 0. 



1 
1 



(B2) 
(B3) 
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We assume that the density matrix giving the initial state induces no spatial variation in X 
at later times (see examples in jisj). In addition, the tree-level kinetic equation Eq. f IBSp 
implies that interactions do not induce further X dependence dxG until at least order ecou, 
so in the constraint equation Eq. ( 1B2I) we can drop dxG as it is O^el^n). These equations 
of motion and translational invariance imply that the most general solution for the Green's 
functions for if at tree level must be of the form [i^ (cf. Eq. (133|) ) 

G^>{k;X) = 2Tr6{e - m^)[e{k''){l + f{k)) + e{-k'')f{-k)] (B4a) 
G'°<(A:; X) = 27r6{k^ - m')[e{k^)f{k) + e{-k^){l + f{-k))\ (B4b) 

for G^, and 

for the retarded and advanced propagators. Note that are proportional to the spectral 
function p = G^ — G^. The functions /, / are required to appear in both Eqs. f lB4al) and 
( ]B4b[) since p = G^ — G^ also. The time- and anti-time-ordered propagators can be found 
by using Eq. ( ]A14I) . _ 

The functions /, / give the distribution of particle and antiparticle states in the plasma. 
Since the right-hand side of the kinetic equation is zero, there is no relaxation or dissipation 
of /, /, which remain static. 



2. First order: Shifted Spectrum 

At first order in y, the self-energy comes from the tadpole diagram with the same topology 
as in Fig. [T^, and is given by 

^^^4/ ^,fG'^{z,w) 



2 V G\{z,w)J 

or in momentum space, 

which is diagonal since / (i^pG^(p) = / d^pG^j^{p), and in fact is independent of k. The 
propagator for A is given by the form Eq. ( 1B4|) with a distribution function Ja = J'a- We 
then insert this into the right-hand sides of the constraint and kinetic equations Eqs. ( lAlSaP 



and (lAT3bl) . Since we can use the tree-level expressions for all Green's functions on the 
right-hand sides, the O operators give zero when acting on them. Furthermore, the right- 
hand side of the kinetic equation vanishes identically since Hi is diagonal (in OTP space) 
and commutes with G. The only effect of Hi is to add a correction to the mass term in the 
constraint equation: 

_ ^2 _ ^n^(^)) X) = 2 Q (B8) 

k ■ dxG{k; X) = . (B9) 
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The vacuum part of the mass correction Hi is absorbed into renormahzation coefficients, and 
the /a dependent part gives a medium-dependent contribution to the mass, as in Eq. (!52|) . 
The form of these equations is still the same as the tree-level equations, with the Green's 
functions translation invariant up to 0{y), so the solution takes the same form as the tree- 
level solution, with a corrected mass: 

— m^ljAl — i€k^ 

G'>ik; X) = 2Tx5{e - m'[fA]mk'){l + f{k)) + e{-k')f{-k)] . (BlOb) 
Because of Eq. ( 1B9I) . there are still no collisions affecting f{k; X) to this order. 



3. Second order: Collisions 

At second order in y, the self-energy receives contributions from the graphs with the same 
topology as in Fig. [T|d. 
In position space, 

ft .^_y^(G\{z,wYG\z,w) -G<{z,wYG>{z,w)\ 

2 \G>{z,wfG>{z,w) -G\{z,wfG\z,w) J " ^ ' 

II2 now generates a nonzero contribution to the right-hand side not only of the constraint 
equation Eq. (lAlSap but also the kinetic equation Eq. f lAlSbp . We again drop the diamond 
operators and the dx in the constraint equation, because dx acting on G is already at least 
order y^. 

We can also rearrange the OTP components of the constraint equation into a more con- 
venient equation for the retarded Green's function Gr = G* — G^ (equivalently for the 
advanced Green's function Ga = G* — G^). We then obtain for the constraint and kinetic 
equations, 

[k^ - - iUiiX) - m2R{k; X)]GR{k; X) = i (B12) 
2k ■ dxG{k; X) = [U>{k; X)G<{k; X) - U<{k; X)G>{k; X)] Q "J^ . (B13) 

The poles of the retarded propagator given by Eq. (lB12p give the medium-dependent masses 
and widths of quasiparticle excitations. 

Note that the kinetic equation Eq. (IBlSp implies for Gr, 

k-dxGR{k;X)=0, (B14) 

so the k ■ dx derivative of Gr (and Ga, and therefore the spectral function p = Gr — Ga) is 
at least third order in y. This will allow us to neglect k ■ dx acting on the spectral function 
in solving the kinetic equation. 

Since the form of the 0{y^) equations of motion Eqs. flB12p and flBlSp is no longer the 
same as the tree-level form, we must be more careful in justifying the form Eq. (IBlOp for 
the Green's functions. The solution to Eq. ( ]B12p for the retarded propagator. 



k^ — — iUi{X) — iU2R[k] X) 
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and similarly for the advanced, gives the spectral function p(fc; X) = {Gr — GA){k; X). The 
spectral function also satisfies p = — G^. Thus, the most general form for G^ is 

G>{k- X) = p{k; X)[l + F{k; X)] + c{k; X) (B16a) 
G<{k; X) = p{k; X)F{k; X) + c{k; X) . (B16b) 

Since the Green's functions satisfied the form Eq. fIBlOl) at 0{y), the function c{k] X) must 
be explicitly at least of order y"^. It must also satisfy the constraint equation. Since c already 
has an explicit y"^ it must satisfy the tree-level constraint equation. Then it is proportional 
to (5(/c^ — m^), and can be absorbed into the p{k;X) part of Eq. ( ]B16[) . Redefining F, we 
bring G^ into similar form as Eq. flBlOp in terms of /, 

G>ik; X) = p{k- x)[e{k'){i + f{k- X)) - e{-e)f{-k- X)] . (bi?) 

Now we evaluate the kinetic equation Eq. flB13p using the expression Eq. fIBlip for 112 
and the form Eq. (IBlOp for the Green's functions. Integrating over positive k^, we find the 
Boltzmann equation for the particle distribution function, 

2k ■ dxfik; X)= — [Cann(A;; X) + C,catt(fc; X)] , (B18) 



where the annihilation and scattering collision terms are the one-flavor analogs of Eqs. ( I59|l 
and dsn]), 

X {f{k)f{k')[l + fA{pW + fA{p)] - [1 + + f{k')]fA{p)fA{p')] , 

(B19) 



and 



oo 



dk'^ n n v\ f d^k' d^p d^p' /n_\4c-4 



Zir zojk J (27r)'^2(X'k' [ZTryzep {Z-Kj-^zep' 

X {f{k)[l + fik')]fAm + fAip')] - [1 + fik)]fik')[l + fAip)]fAip')} . 

(B20) 

We have performed the frequency integrals in the collision term using the free spectral 
function; corrections to this approximation are 0{y^). 

Note that Eq. (IBISP tells us that collisions create 0{y'^) spacetime variations in the 
distribution functions. So if we were to continue to higher orders in perturbation theory, 
we could not keep neglecting the d\ term in the constraint equation Eq. ( lAlSal) and the O 
operators in the constraint and kinetic equations Eqs. (lA13aP and (lAlSbp as we have done 
so far. 



4. Equilibrium Solution 

To find the equilibrium solution for the distribution functions for ip and A, we require 
that the collision terms in Eqs. flBlOp and flB20p vanish. It is easy to show that, due to the 
property 

1 + nB{ko -p)= e^'^"-'^)/^ - p) , (B21) 
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and the momentum-conserving delta functions in Eq. (IBlSp . the forms 

/(fc;X)=n^(A;o-/i(X)), /(A;; X) = n^(A;° + /i(X)) 

fAik;X)=nBiko), ^ ' 

are such a solution, as long as the X variation of the chemical potentials is smaller than 
0{y'^). Therefore, interactions at 0{y'^) (namely, binary collisions or annihilation due to 
interactions) leave the distributions Eq. (]B22[) unchanged, as long as satisfies the 

continuity equation k ■ dxfJ^ = 0. 

Appendix C: Two-Flavor Spectrum of Excitations to 0{e) 

In this appendix we include further details on the two-fiavor CTP Green's functions and 
solve for the spectrum of excitations to 0{e) in the toy model introduced in Sec. [Tll 

The spectrum of excitations can be found from the equations for the Wigner-transformed 
retarded and advanced propagators, 

- G^'^{k- 1) + \e-^^{m\t),G^'^{k- 1)] - ze-^^{A:°S(t), t)] ^^^^ 

and 

2fc°aiG^'^(fc; t) + le-'"" [m\t), G'^'^ik; t)] + ^e"'* [A;°S(t), t)] = e"'^ [H^, . 

(C2) 

After solving for G^'^ and constructing the spectral function p = G^ — G^, we can then use 
p = G^ — G^ and the kinetic equations Eq. (l26l) for G^ to determine G^. 
We begin with the constraint equations at 0{e^): 



k'G--[m\G^=z r , (C3) 



proportional to the identity in both field space and CTP space. We gave the general form of 
the solution for Gfj in Eq. (133|) . The solution for the diagonal components of the retarded 
and advanced propagators are 

^11% = 75 2 I • ; ■ (^^) 

k-^ — mi 2 ± i^k^ 

Since all CTP components of the off-diagonal Green's function G12 obey the same equation 
Eq. (ICSp . and since G^'^ = G* — G^, this implies that the 0{e^) off-diagonal retarded and 
advanced propagators Gi2'^ = 0. 

Next, keeping terms up to 0{e) in Eq. fICip . we obtain the equation for retarded and 
advanced propagators, 

^2^H,A _ 1 1^2^ ^R,Aj + .k'^iEit), G«'^} = Z + ^{n^'^, G^'^} . (C5) 
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Writing the diagonal and off-diagonal components (in field space) separately, we find 

{k^ -ml + 2ik^T.u - 2infi)Gfi = i + iili^^ - A;°Si2)Gfi + (n^^ - k^^2i)G^2 (C6a) 



- + ml) + ^^°(Sn + S22) - ^(nfi + nf^) ]G^2 = ^(nfs - A;°Ei2)(Gfi + G^2) 



2 



(C6b) 



and similarly for G22 and (and G"^). 

In the equations for the diagonal components, the off-diagonal components G12 21 si-ct as 
sources for G'ii,22- As they multiply factors 11, E explicitly of C(e), we can use the C(e°) 
solutions for G{^2 21 these terms, which we noted above are zero. Therefore, the 0{e) 



solutions for Gi{\2 simply 

~ F - 2(t) + 22A:0Sn,22(t) - 2^nfi';^2(^; ^) ' 

These solutions give us the shifted spectra of the excitations that had tree-level masses mi^2 
due to interactions S with the wall and 11 with other particles in the bath. They induce 
a corrected mass and width for the excitations. From these solutions, first, we construct 
the diagonal spectral functions pii,22 = Cf^ 22 ~ 22' ^^^^ second, we can use the spectral 
functions to construct the Green's functions G\ 22- 

G> ^22(^; t) = pn,22(fc; t)[^(^°)(l + /ii,22(A;; t)) - 9i-k^)fnM-k; t)] (C8a) 
Gfi,22ik; t) = puMk; t)mh,Mk; t) - e(-fcO)(l + /ii,22(-A;; t))] . (C8b) 

These Green's functions have the same form as the tree-level solutions Eq. fl5^ . with mod- 
ified dispersion relations given by the spectral function. These forms satisfy the constraint 
equations Eq. fl25l) for Gf^^ 22 ^i^)^ while the kinetic equations Eq. fl26l) at 0(e) give non- 
trivial evolution of the distribution functions /, /, as we derived in the main text. In this 
paper, since we were interested in the evolution of /, / only at 0{e), we kept the tree-level 
spectral functions given in Eq. fl33l) when solving the kinetic equation. 

For the off-diagonal components, substituting the solutions Eq. (1C7P for Gf{'^2 i'^to the 
right-hand side of Eq. (lC6bp implies that, to 0{e), 

r^,A.,.._ 2[nf^-^(fc;t)-fc°Si2(t)] 

(P-m?(t))(P-m2(t))- ^^^^ 

This solution induces an C(e) contribution to the spectral function pi2 = Gf2 — G12, which 
would add a modification to the dispersion relations appearing in the tree-level solutions 
Eq. f l55]l for Gf2, so that these also obey the constraint equations Eq. fl2S]) at C(e). However, 
to solve for the C(e) evolution of /i2,/i2 with the kinetic equation, it suffices to keep the 
tree- level dispersion relations in Eq. ( 133|) . 
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